


nterLINES 





Thermal rub effect in rotating machines 





by Agnes Muszynska, PhD 
Senior Research Sctentist 
Bently Rotor Dynamics Research Corporation 


rotor-to-stator rubbing in heavy rotating machines, A 
simple model explains the often-observed-in-the-field 
“spirals” in the 1X response at a constant rotative speed. 


T his article discusses the thermal;mechanical effect of 


Introduction 


Rotor-to-stator rubbing in rotating machines is a malfunc- 
lion associated with a physical contact between rotating and 
stationary parts, which, during normal operation, should not 
be touching, This contact causes changes in the system force 
balance and dynamic stiffness and results in modifications of 
the machine motion. The effect is usually associated with a 
decrease in the energy provided to maintain the main opera- 
tional motion (rotational mode}, and an increase in the level 
of “parasitic” vibrations, 

A rolor-to-stator rub is a secondary malfunction resulting 
from a primary cause which perturbs the machine during 
normal operational conditions. This primary cause can origi- 
nate from various sources, such as unbalance, thermal and/or 
assembly misalignments, stator/casing motion, fluid-induced 
rotor self-excited vibrations due to fluid forces in bearings, 
and/or in the main flow, and in seals, in particular. The 
occurrence of the primary source is manifested by changes in 
rotor centerline position and/or shaft lateral vibrations so that 
allowable clearances are exceeded, and the rotor-to-stator 
contact occurs. 

A rotor-to-stator rub involves several physical phenomena, 
the most important being friction, impacting, coupling and 
stiffening effects, 
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The rub-related rotor-to-stator friction is responsible for 
the rotor thermal balance change, known as “Newkirk effect.” 
Based on Taylor’s results [1], Newkirk [2] pointed out that, 
when a rubbing rotor is running below its first balance reso- 
nance speed, the rub-induced rotor lateral vibrations tend to 
increase in time. Later, this effect was studied by many other 
researchers, who confirmed that these vibrations can grow in 
amplitude and phase, resulting in “spiral” vibrations. 


Thermal effect of friction 


During machine operation, the rub of an unbalanced shaft 
usually takes place at the seal with the smallest clearance 
which is closest to the shaft antinodal position. Rotors can 
also rub against packings or oil deflectors. The rub occurs at 
the shaft “high spot” (angular location under the highest 
iension stress). At a constant rotative speed, for the mode of 
vibration predominantly synchronous with a circular or 
slightly elliptical Orbit, the rubbing high spot occurs at the 
same shaft location at each turn. The rub at the high spot 
causes friction-related heating and local thermal expansion 
of the shaft. Due to normal or accidental presence of fluids in 
the rotor/stator clearance areas, the rub-generated heat can 
be carried away by the fluid flow, so that the shaft heating and 
thermal expansion process may be relatively slow. Due to 
local expansion, the shaft bows, causing an additional unbal- 
ance in the rotor, 

At a rotative speed lower than the first balance resonance, 
the phase angle between the high spot and the heavy spot 
(angular location of the unbalance) is less than 90 degrees 
(Figure 1); the additional bow-related unbalance then adds to 
the originally existing one (vectorially), and causes an 
increase of the overall unbalance inertia force acting on the 
rotor. The relationship between the new unbalance and the 
original one explains the phase lag growth, as the phase 
between the force and response adjusts to the value deter- 
mined by the system parameters and a constant rotative 
speed. The phase lag increase causes the high spot to change, 
thus the rubbing location changes around the rotor circum- 
ference. The increase of total value of unbalance causes 
higher vibration amplitudes, which may result in stronger 
rubbing and larger amounts of heat. LE the rub-related heating 
is fast enough, the rotor local thermal expansion becomes 
high and may exceed the elastic limit. The rotor then becomes 
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permanently bowed due to local plastic deformation and can 
no longer be operated. 

More often, however, the process of rub-related, local 
thermal expansion and shaft bowing is slow, allowing for 
continuous phase adjustments. The resulting synchronous 
vibration amplitude changes are also slow The phase lag 
growth may exhibit several full circles. This increasing phase 
lag reflects the continuous change of the high spot rubbing 
location. According to field reports, the periods of full circle 
phase lag growth may vary from a few minutes to over 24 hours 
(Figure 2). 


Simple mathematical model of the thermal rubbing 
phenomenon 


The effect of the rotor-to-stator thermal rub that is 
observed most often in rotating machines is the slowly 
increasing (or fluctuating) 1X vibration amplitude, with a 
definite growth of the 1X phase lag (Figure 2). The effect of 
this 1X amplitude variation and the phase growth can be 
explained using a simple mathematical model. 

Consider the isotropic (symmetric) rotor of a fluid- 
handling machine at its first lateral mode. The rotor is rub- 
bing against the stationary part at an axial location “R” The 
rubbing-related, local thermal expansion causes the rotor to 
bow, This results in an additional unbalance. The rotor model 
is as follows: (1) 


M2 + DZ + Kz+ DG janz) =mrořel * 8) 4 xkp, (pelle as! 


d 


a =arctan (y/x) Z=X+f¥y j= = 


where z is the rotor lateral displacement combined in 
ane complex number, x is horizontal, y 1s verti- 
cal displacement of the rotor in stationary coor- 
dinates 

KMO are rotor modal stiffness and mass, respectively 
D. is external damping 


D is the rotor surrounding fluid radial damping 

A is the fluid circumferential average velocity 
ratio (lambda) 

n is rotor rdtative speed (omega) 

m., are the rotor modal unbalance, mass radius and 
angular orientation (delta), respectively 

p_(t) is Tub-related, thermally-induced shaft bow 
expressed in coordinates rotating with the 
shaft. This bow is a function of time t (rho sub -) 

K is a modal factor which depends on the axial 
location of the rubbing area (kappa) 

a Isa constant angle related to the shaft twist 


between the rubbing section and the modal 
mass location (alpha sub s) 
The thermal bow occurs at the angular location of the high 
spot (with a small shift due to friction angle). The new ther- 
mal, bow-related unbalance adds vectorially to the original 
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Figure 1 

At a constant rotative speed, 2, below the first balance resonance, 
the phase lag, ay, between the heavy spot (original unbalance 
force) and high spot is less than 90° (a). Due to the rub, the shaft 
bows in the direction of the high spot and a new unbalance force, F 
occurs {b}. The original and rub-related unbalances add together, 
producing the “effective” unbalance force {c}. At a constant rota- 
tive speed, %, the phase lag between the effective force and 
response (high spot) must be constant. This means that the shaft 
has to rotate, yielding a new high spot (d). Note, also, the increase 
of Orbit magnitude at the end of this process. 


i el a ery 










vn ETA EEE i SIENET 


"P 
WALI: bat Jrg: gi TF " 


ainia y aL Bad i es del ges 
Pas E Tharp eT] BLEJ? To w hI BELT 
ei all = 
Q MI ' ; 
de “gai dey aa PTE thal 
a Eel a Osir EL | JLAM EJ 
aJ bEJILT aeg i i 
s ETETE i : a a Os de 
=p Os = ee J: aaarit 
PETE, a Gon EET a = -ag Iak 7 FALDER wa 
i ; A manir Prat ajiak damea Wi 
| i ae Nee lee e aa | aaenttay 
i s afar Da naa PIDE FT re 
Srdbie ty as coy ge | ee 


LERE] 
CRER HI 


i 
| jinzi eee 
| esi bmi jMardbsdegardarda gy 















y EL an [a 
Paaa Jae Ta 
alae d a: 
: aadsai et 


A : l 
l LP H idT 4. | . i i 
-E5 I $ 4J gsl Jz 

i tI: Eri AT beau aa: #4 ae = . baa y sas as Í i 

i Seon a] | 

eh, eds bay Eee Tapie i et Hale bi 

| | RANE i eee " m paap ia debe i ao) Garay 
i vizet EE | AJILETA aace s "a a i 


í 4 aga setae 2 peal dita ae i 7 afi aint ica die : aeia 
ko amena L S See EEL f F 
i aa War aT amnaid 
i LEHE HEF ; i É, 
:G4: a . CHG ee ki 
w EEN IR biiga] 


#3: 6: 1 LEELEE a 
saraa a — enor ist yy 
“#2 Bearing ¥ Probe X° 


ee ae 
Steady State 1X Filtered Uncomp 


Rub occured al an oil seal in the 
Thrust Bearing Box near Bearing #1 


3.0 mil pp Full Scale CW Rotation 


Figure 2 
Thermal bow-related rub diagnosis of a turbine rotor: Polar plot of 
time-trended synchronous (1X) response vector, Rotative speed & 
= 3600 rpm. During 30 minutes, the 1X vector accomplished six 
full circles. 
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unbalance. When the rotor runs at a speed below the first 
balance resonance, the angle between the heavy spot (origi- 
nal unbalance) and high spot (direction of the new rub- 
related unbalance) is less than 90 degrees. The vectorial sum 
of them is, therefore, higher than the original unbalance, and 
the resulting phase angle has to adjust to the phase appropri- 
ate for the specific constant speed of rotation. Figure 1 illus- 
trates this effect. The model Equation (1) explains the 
phenomenon mathematically, 


Assume that the thermal bow p-p is a known, slowly varying 
function of time. Assume, also, that there was no rubbing 
during startup of the rotor up to the rotative speed R. At the 
instant just before the rub occurred, the rotor vibrations 
described by Equation (1) were synchronous due to original 
unbalance only {or (0) =0): 





z(t) | p= Aget 6 + aa) (2) 
where the amplitude and the phase are as follows: (3) 
mra” lD + D{(1-A)] 
Å= — ) 8g =arcan———____— 





Mark 


(K-Mo}? + 92(D, + D(1-A)]? 


The angle 4 + œp deseribes the high spot orientation, as 
measured by the Keyphasor™ transducer. 

At the moment when the rub occurs, the synchronous 
vibrations start changing because the force balance, Equation 
(1). acquires a new, synchronously rotating force related to the 
thermally-induced shaft bow unbalance: 


F-K p_(thell™ to-as+aq)  becausea=Qt+ő+ap 


Z|y= Agel OF 40) 4 (A, j jett iasta (4) 
= g elSt+ & + agt) 
where 


«K(or j} 
AY |; = 





(K-Ma2)2 + o2 [Da + D(1-a)]2 


is the amplitude of the response to the shaft bow-related 
unbalance force. The corresponding response phase is the 
same, cy, 48 defined previously. Two elements of the response 
(4) can be vectorially added, The newly-acquired amplitude, 
B, and phase 8), of the rotor total synchronous response are 
as follows: 








By=V AZ+IAS] 2+ 2A | JAgcos(ag—as) (5) 
tan(a,-cs) 
G4 = arctan — 
Ag 
1+ 


(A; ,}cos(ap-ag) 
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Figure 3. 
Schematic of a rotor rubbing at the axial location, R. 


It can be seen that, in comparison to that of Equation (2), 
the synchronous amplitude B, in Equation (4) is magnified by 
some factor, and the phase lag becomes higher by the amount 
Pr 

At the next moment of the continuously rubbing rotor, 
{these moments, denoted by subscripts | ,,i=0,2,..., belong 
to a different, much slower, time scale than t) the synchronous 





vibrations are: (6) 
x(t) >= Agl to +0) 4 (a, jellet+ dag + ety + 84) 
i a 
= ByellAt+ 5 +19 + Bo) 
where 
(erit |a) 7 
Ajaj — 2 (7) 
(erit |+) 


and the newly-acquired LX amplitude and phase are: 








Ba= y AR+(A) + 2A (Ay | Jopas — (8) 
tanigani + 8 
È» = arctan a (gas + 81) (9) 
Ag 


(A, | JCOS(ag-ag + By) 

At this moment, further amplitude and phase lag growths 
are noted, 

By continuing the calculating process, at any next moment 
“i the rotor synchronous vibrations can be calculated based 
on the previous moment “i-l; and will have the following 
form: (10) 

j(at+ 6 +a) jiat+ -0g + 2a_,+ Bi 
z(t) |= Ago! 0) + (a Jolt ras + 2a0 + Bia} 


= Biellatt 6 + Cig + Bi), = joes ee 
where 


| Kletti] 
Ajj= - (11) 


fiK-Min2)2 + 92[D. + D(1-A)]? 











AS + (Ay | JP-+ 2Ap(A, | JOOSlag-as + Bj._4 
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tanlag-as + j4) : 
Bi =arctan se (13) 
(A, | Jcoslag-ag + 84) 

Equation (10) describes the rotor spiraling/oscillating 
mode of synchronous vibrations, The full rotation of the 
phase is reached for the value of “1” when §; = 360°, The 
number of steps “i” in the slow time scale needed to accom- 
plish the full rotation depends on the tnitial phase lag, ap, 
which, in turn, depends on the closeness of the rotative speed 
to the first balance resonance frequency, YKM. The closer the 
valuc of @ is to YKM, the higher the original phase lag ag is, 
and the faster the full rotation is accomplished. Numerical 
examples below illustrate this effect. 

Numerical examples 

Example I 

In this numerical example, it is assumed that the original 
unbalance-to-thermally-induced unbalance amplitude ratio 
Is constant: 


Ky 
Assume that a, = —40°,6 = 0, a, = 1°. 
The rotative speed, Q, to the system natural frequency ratio 
can be expressed in terms of the phase lag am and the damp- 
ing factor g = [De + D(1-A)|/2VKM | 
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0.02, the frequency ratio is 














For the damping factor Z = 
WVEM = 0.98. 

Calculation results are presented in Table 1 and Figure 4. 

As can be seen, the full 360° rotation of the response phase 
was accomplished in twelve steps (1). After an initial increase 
between steps O and 1, the response amplitude 1s slowly 
variable, 







a rT 


Oy + È; 





=-141497 +V 1+ (1.192) 


ies e e a pe epa 





-270° 


Ait, = 24 
_ alpha, = 1° 







alpha, = -40° 






-180° 


-90" 
Figure 4 
Polar plot of the rotor synchronous (1%) response, One full rotation 
of the phase accomplished in twelve steps “1” {Example 1). 
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Figure 5 
Polar plot of the rotor synchronous (1X) response, One full rotation 
of the phase accomplished in six steps “i” (Example 2). 












Table | 


Rubbing rotor to non-rubbing rotor 1X response amplitude ratio and its phase as functions of “slow” time steps “T” 
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Example 2 For the damping ratio g = 0.02, YEM = 0.99. 

In comparison to Example 1, in this example the phase angle 
lag, &p 18 higher. That is, itis assumed that the rotative speed is 
closer to the resonance. 


‘The results are presented in Table 2 and Figure 5, 


In this example, the full 360° rotation of the response phase 
was accomplished in less than six steps (1). The Polar plots for 


tig = 607 both Examples 1 and 2 are very similar; they differ, however, by 

F TETE the time period of the full rotation. The closer the rotative 

— = 0.587 + V1 + (0.58g)? speed is to resonance, the faster the full phase rotation is 
VKM accomplished. 
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Table 2 
Rubbing rotor to non-rubbing rotor 1X response amplitude ratio and its phase as functions of “slow” time steps “i” 


Example 3 
In this example, the rub-related thermal how unbalance is 


aol on o ; F The remaining parameters are the same as in Example 2. 
considered increasing, according to the following relationship: | 


The results are presented in Table 3 and Figure 6. 





kKpy] | The full 360° phase rotation is accomplished in less than six 
=A fA ,=2+011 steps, similar to Example 2. The amplitude is slowly increasing, 
mro demonstrating the spiraling effect. 
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Table 3 
Rubbing rotor to non-rubbing rotor 1X response amplitude ratio and its phase as functions of “slow” time steps “i” 
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Figure 6 
Polar plot of the rotor synchronous (14) response for the case when 
the rub-induced bow is an increasing function of time (Example 3). 


Conclusions 


The model of the thermal rub phenomenon presented in this 
article is certainly very simplified, It was assumed, in particular, 
that the rub occurs al some axial location, and its effect on the 
first mode symmetric rotor behavior was limited to rator bow- 
ing. The amount of the rotor thermally-induced bow was not 
quantified, nor correlated with the rubbing rotor-to-stator con- 
tact frictional and local expansion effects. [1 was only quali- 
lalively shown that, at a constant rotative speed, the thermally- 
induced rotor bow, as a result of rubbing, creates additional 
unbalance, and, consequently, causes phase lagging as the 
rubbing spot continuously changes. 


The assumption about the rotor symmetry (isotropy) and 
limitation to only one lateral mode provides another simplifica- 
tion. Both of these latter assumptions can easily be alleviated 
without losing any qualitative features. In particular, if the 
rotor speed is in the range of the second mode, slightly lower 
than the second balance resonance, and if the rub occurs at the 
rotor axial location in the area of the highest unbalance, then 
the angular situation between the original unbalance and rub- 
related unbalance is practically the same as in the simple one- 
mode example considered above. Therefore, the rubbing rotor 
behavior will be very similar, as discussed in this article. W 
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Notation 
Rotor synchronous response amplitude with and 
without thermal rub effects 
Rotor synchronous response amplitude affected by 
rub thermal effect 
Fluid damping and shaft modal external damping, 
respectively 
Integer describing time instants 


Rotor modal stiffness 

Unbalance, mass radius and angular orientation, 
respectively 

Kotor modal mass 

Time 


mass location 

Rotor lateral coordinates at the modal mass 

Rotor synchronous response at instant “i” 

Rotor synchronous response phase without ther- 
mal effect 

Rotor instantaneous angular orientation at modal 
mass 

Angular twist between the rubbing place and 
modal mass location 

Modal factor 

Fluid circumferential average velocity ratio 
Rub-related thermally induced shaft bow 
expressed in rotating coordinates 

Rolative speed 
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